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ORIGINAL ARTICLE 

Oncogenic B-Raf signaling in melanoma cells controls a network 
of microRNAs with combinatorial functions 

KL Couts 1 , EM Anderson 2 , MM Gross 2 , K Sullivan 2 and NG Ann 1 ' 3 

Over two-thirds of melanomas have activating mutations in B-Raf, leading to constitutive activation of the B-Raf/MKK/ERK signaling 
pathway. The most prevalent mutation, B-RafV600E, promotes cancer cell behavior through mechanisms that are still incompletely 
defined. Here, we used a sensitive microarray profiling platform to compare microRNA (miRNA) expression levels between primary 
melanocytes and B-RafV600E-positive melanoma cell lines, and between melanoma cells treated in the presence and absence of an 
MKK1/2 inhibitor. We identified a network of >20 miRNAs deregulated by B-Raf/MKK/ERK in melanoma cells, the majority of which 
modulate the expression of key cancer regulatory genes and functions. Importantly, miRNAs within the network converge on 
protein regulation and cancer phenotypes, suggesting that these miRNAs might function combinatorially. We show that miRNAs 
augment effects on protein repression and cell invasion when co-expressed, and gene-specific latency and interference effects 
between miRNAs were also observed. Thus, B-Raf/MKK/ERK controls key aspects of cancer cell behavior and gene expression 
by modulating a network of miRNAs with cross-regulatory functions. The findings highlight the potential for complex interactions 
between coordinately regulated miRNAs within a network. 
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INTRODUCTION 

Oncogenic mutations in the protein kinase, B-Raf, occur in 50-70% 
of sporadic melanomas. 1,2 The most frequent mutation, 
B-RafV600E, elevates kinase activity and constitutively activates 
the MKK/ERK signaling pathway, which promotes tumor 
progression and metastasis through enhanced cell proliferation, 
survival, motility and invasion. 3,4 Although the prevalence of the 
B-RafV600E mutation and cellular responses have been intensively 
studied, molecular targets of this pathway in melanoma are still 
incompletely understood. 

MicroRNAs (miRNAs) are short, non-coding RNAs that provide 
feedback mechanisms to maintain stable protein expression and 
confer robustness in development and normal cell homeostasis. 5 
miRNAs also modulate protein levels and cell phenotypes in 
response to environmental stress, and are often mutated or 
deregulated in human diseases. Widespread miRNA deregulation 
has been reported in cancers, 6 and certain miRNAs have been 
classified as oncogenes or tumor suppressors across multiple 
cancer types. Although changes in miRNA levels are not typically 
sufficient for promoting tumorigenesis, miRNAs are necessary for 
controlling cancer cell properties and maintaining disease. 

Several oncogenic signaling pathways have been shown to 
support tumorigenesis by modulating miRNAs. Altered transcrip- 
tion factor expression is one mechanism that has been linked to 
the deregulation of cancer-associated miRNAs. For example, c-Myc 
induces the oncogenic miR-17~92 cluster and represses other 
miRNAs, including let-7, miR-22, miR-26a and miR-34a, in B-cell 
lymphomas. 10 In neuroblastomas, the related transcription 
factor, MYCN, promotes oncogenesis by inducing miR-17~92 



and miR-221. 11 Two studies have investigated the correlation 
between B-Raf mutational status and miRNA expression in 
melanomas, and only one was able to link three miRNAs to 
B-RafV600E. 12,13 Other studies have reported large numbers of 
differentially expressed miRNAs in melanomas vs normal 
melanocytes, but little agreement was observed between 
them. 14-16 As a result, the miRNAs that are deregulated in 
melanomas, and specific targets of oncogenic B-Raf, have yet to 
be established. 

In this study, we used a highly sensitive microarray profiling 
platform to identify miRNAs regulated by oncogenic B-Raf 
signaling in melanoma cells. We show that a network of >20 
miRNAs are controlled by B-Raf/MKK/ERK signaling, encompassing 
more than half of the miRNAs which are deregulated in melanoma 
cells compared with primary melanocytes. Cell treatment with 
miRNA mimics and inhibitors reveals that the majority of these 
miRNAs control cell proliferation and/or invasion. Target predic- 
tion algorithms suggested that many cancer-associated genes are 
targets for convergence by miRNAs within the network. This was 
validated at the level of protein expression and 3'UTR activity, 
raising the possibility that miRNAs functionally interact with one 
another. Both cell invasion and protein repression were augmen- 
ted upon combinatorial miRNA expression, and more complex 
interactions between miRNAs, such as interference and latency 
effects, were observed for protein expression. Altogether, 
our findings establish a network of miRNAs under the control 
of oncogenic B-Raf in melanoma and highlight the complexity 
of cross-regulatory interactions between miRNAs in modulating 
cancer cell properties. 
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RESULTS 

MicroRNA profiling in primary melanocytes and established 
melanoma cell lines 

Microarray profiling was used to identify miRNAs both deregu- 
lated in melanoma cells and responsive to signaling through 
B-Raf/MKK/ERK. We analyzed miRNA expression levels in two 
primary melanocyte cell lines (NHEM2493 and HEM-LP) and 
compared these with six melanoma cell lines established 
from human primary vertical growth phase tumors (WM115, 
WM278 and WM793) or metastatic tumors (WM239A, 
WM1617 and 1205Lu). To distinguish miRNAs responsive to 
constitutive B-Raf/MKK/ERK signaling, we carried out an indepen- 
dent experiment comparing the six melanoma lines treated 
with and without the MKK1/2 inhibitor, U0126 (Figure 1a). 
All melanoma cell lines carried activating B-Raf-V600E or 
-V600D mutations 1,1 7,1 8 and expressed elevated levels of 
phospho-ERK1/2, which was attenuated by U0126 treatment 
(Supplementary Figure S1). 

Biological replicates were analyzed using microarrays contain- 
ing triplicate probe sets of 470 miRNAs (Thermo Fisher Scientific, 
Lafayette, CO, USA; miRBase 9.0). An miRNA reference probe set 
was added to each sample, allowing better comparison of features 
within a single array and across many arrays. In total, 18 800 
measurements were generated and normalized for variability in 
labeling efficiency, hybridization and baseline fluorescence 
(Supplementary Table 1). Of the miRNAs surveyed, 171 were 
excluded from further analysis based on their score of 'present' 
(P<0.05) in fewer than 10 of the 40 samples (Supplementary 
Table 2). 

Averaged intensity measurements identified >30 miRNAs 
which differed between melanocyte and melanoma cell lines by 
1.5-fold or more (P<0.05, Figure 1b), as well as an additional 
24 miRNAs which differed by at least 1.3-fold (P<0.05, 
Supplementary Figure S2). Quantification of several miRNAs was 
confirmed by real-time qRT-PCR. The results showed excellent 
agreement with microarray experiments across all cell lines 
(Figures 2a,c), validating the magnitude and directionality of 
deregulation for each miRNA. 



B-Raf/MKK/ERK regulates a network of cancer-associated miRNAs 
To distinguish which miRNAs are targets of B-Raf/MKK/ERK 
signaling, we examined changes in miRNA expression following 
treatment with MKK1/2 inhibitor (Figure 1b). More than half of 
miRNAs altered by ^ 1.5-fold in melanoma cells vs melanocytes 
also showed consistent, reciprocal changes in expression after 
U0126 treatment. Several family and cluster members of these 
miRNAs were also regulated by U0126 but were only weakly 
regulated (1.3-fold or less) in melanoma vs melanocyte cells 
(Supplementary Figure S2). Although U0126 elicited relatively 
small changes in miRNA expression, the measurements were 
validated by qRT-PCR analysis of cells treated with U0126 as well 
as CI-1040, another MKK1/2 inhibitor which also decreases 
phospho-ERK1/2 levels in melanoma cells (Supplementary Figure 
Sib). Both U0126 and CI-1040 demonstrated consistent changes 
in miRNA expression compared with microarray measurements 
(Figures 2b and d). A previous study in HeLa cells showed that 
miR-17 and miR-20a were increased, and let-7 decreased each by 
20% in response to U0126, 19 consistent with their responses to 
drug in our study. 

Most B-Raf-responsive miRNAs identified in our screen have not 
been extensively studied in melanoma, but show oncogenic or 
tumor suppressor roles in other cancers (Table 1). Among the best 
studied is miR-17~92, which is elevated in many cancers and 
associated with increased proliferation and tumorigenesis. 20 
Conversely, miR-34a and let-7 family members are tumor 
suppressors in cancers, and inhibit proliferation by modulating 
key cell-cycle regulatory genes. 21,22 Together, the results reveal a 
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Figure 1. B-Raf/MKK/ERK regulates miRNAs differentially expressed 
between melanocyte and melanoma cells, (a) Strategy for identify- 
ing B-Raf/MKK/ERK-regulated miRNAs in melanoma cells, (b) Heat 
map indicating fold changes in miRNA expression between 
melanocytes ('Mcytes') vs melanomas, and dimethyl sulfoxide 
(DMSO)- vs 110126-treated cells. Signal intensities averaged between 
biological replicates for each melanoma cell line were normalized to 
the averaged signal between two melanocyte cell lines. Averaged 
signal intensities of 110126-treated cells were normalized to the 
average intensity of the corresponding DMSO-treated cells. Data are 
represented on a log 2 scale. 

network of cancer-associated miRNAs that are controlled by 
oncogenic B-Raf signaling in melanoma. 



B-Raf/MKK/ERK-regulated miRNAs control cell proliferation and 
invasion 

To investigate cellular functions of miRNAs responsive to B-Raf/ 
MKK/ERK, we examined effects on cell proliferation and invasion. 
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B-Raf/MKK/ERK downregulated miRNAs 
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Figure 2. Real-time PCR validation of microarray data. Quantitative RT-PCR for selected miRNAs was performed and compared with 
microarray results, (a, c) miRNA expression across cell lines from microarray and qRT-PCR analyses. Expression in each cell line was normalized 
to the averaged signal from two melanocyte cell lines, (b, d) Percent change in signal after U0126 treatment from microarray and qRT-PCR 
analyses, and CI-1040 treatment from qRT-PCR analyses. Values represent average ±s.e.m. from two combined biological replicates 
(microarray) or three technical replicates (qRT-PCR). The qRT-PCR experiment was performed independently using samples prepared 
separately from the microarray experiment. 



Synthetic oligonucleotide miRNA mimics were transfected into 
WM239A cells, and proliferation was quantified using a colori- 
metric assay. Of the 25 miRNAs tested, 14 altered cell proliferation 
in a manner consistent with their direction of regulation by 
signaling (Figure 3a). Seven miRNAs that were inhibited by B-Raf/ 
MKK/ERK in melanoma cells suppressed proliferation, including 
let-7i, miR-34a and miR-22, which repressed proliferation in other 
cell systems (Table 1). Likewise, seven miRNAs induced by B-Raf/ 
MKK/ERK increased proliferation, including miR-17 and miR-92/25 
family members, which enhanced proliferation and tumorigenicity 
in other cancers (Table 1). Although the magnitudes of 
proliferative responses in the 96-well assays were modest, 
transfection by miRNAs in each case led to qualitative differences 
in cell confluency compared with control treated cells, after 
several days (Figure 3b). 

We next examined the network of miRNAs for effects on cell 
invasion using a three-dimensional spheroid assay. Spheroids 
were formed from melanoma cells transfected with miRNA mimics 
and embedded into collagen, allowing growth and invasion to be 



monitored in an extracellular matrix with composition similar to 
that of human dermis. Striking phenotypes were observed, where 
miRNAs induced by B-Raf/MKK/ERK significantly enhanced inva- 
sion into matrix (Figure 3c; Supplementary Figure S3). On the 
other hand, miRNAs that were downregulated by B-Raf/MKK/ERK 
showed mixed effects. Some inhibited invasion (let-7i, miR-211 
and miR-22; Figure 3d; Supplementary Figure S3), suggesting that 
their repression by B-Raf signaling promotes invasion in melano- 
mas. Others increased the extent of cell invasion (miR-34a and 
miR-26a; Figure 3d; Supplementary Figure S3), suggesting that 
they targeted genes controlling mechanisms inhibiting invasion 
as well as mechanisms promoting cell proliferation. Although let- 
7i, miR-211 and miR-22 decreased cell proliferation, which might 
have accounted for the reduced spheroid invasion, these miRNAs 
also reduced cell motility in transwell migration assays (data not 
shown) suggesting that loss of spheroid invasion was not solely 
due to decreased cell number. Because B-Raf/MKK/ERK often co- 
regulated members of the same miRNA family (for example, let-7i/ 
7g, miR-1 7/20/1 06/93, miR-92/25 and miR-221/222), functional 
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Table 1. miRNAs regulated by B-Raf/MKK/ERK in melanoma cells appear in other cancers 



miRNAs 


Deregulation in cancer* 


Associated functions 3 


Known targets b 


(Downregulated) 




Suppresses proliferation and metastasis 




let-7 


Decreased in lung, breast, prostate, 


CCND1, NRAS, KRAS, HMGA2 




gastric cancer, hepatocelluar carcinoma 


Promotes apoptosis 


CDK6, E2F3, CCNE2/BCL2 33 


miR-34 


Decreased in lung, pancreatic cancer, 


Suppresses proliferation, tumorigenesis 




neuroblastoma 


Promotes apoptosis, senescence 




miR-22 


Decreased in breast cancer 34 


r IT ^- 34 

Suppresses proliferation 


i-r-» 34 

ERor 


miR-26 


Decreased in thyroid cancer 


Suppresses proliferation 


EZH2, PTEN 35 




Increased in gliomas 


Promotes apoptosis 




miR-21 1 


Increased in pancreatic cancer 


Unknown 


Unknown 


miR-125 


Decreased in breast cancer, neuroblastoma 


Suppresses proliferation, 


ERBB2, ERBB3 




Decreased in anaplastic thyroid carcinoma 36 


migration, invasion 


TGFBR1/SMAD2 36 


miR-30 


Suppresses invasion 36 


miR-10 


Increased in breast cancer and glioblastoma 


Promotes invasion, metastasis 


HOXD10, HOXA1 


(Upregulated) 








miR-221/222 


Increased in pancreatic, prostate, cancer, 


Promotes and suppresses proliferation 


KIT, p27, p57 




papillary thyroid carcinoma, glioblastoma 






miR-17-92 


Increased in lung, colon cancer, 


Promotes proliferation, tumorigenesis, 


E2F1, p21, AIB1, PTEN 




B-cell lymphoma, glioblastoma 


angiogenesis 




miR-106a~363 


Increased in T-cell leukemia 


Promotes anchorage-independent growth 37 


p21, RB1, ARID4B 


miR-106b~25 


Increased in prostate cancer, 38 


Promotes proliferation, 


p21, E2F1/BIM 39 




hepatocellular carcinoma 39 


anchorage-independent growth 39 


mTOR 40 


miR-100 


Increased in pancreatic cancer, 


Unknown 




Decreased in ovarian cancer 40 







information from review by Lee and Dutta 6 unless otherwise referenced. b Confirmed targets identified using TarBase database 41 unless otherwise referenced. 



screens used miRNA mimics to avoid complications from 
simultaneously inhibiting multiple family members. However, cells 
transfected with miRNA inhibitors confirmed opposite effects on 
spheroid invasion compared with miRNA mimics (Supplementary 
Figure S4). Together, the results suggest that miRNAs which are 
coordinately upregulated in response to B-Raf/MKK/ERK promote 
the cancer phenotype by enhancing cell proliferation and/or 
invasion, while miRNAs inhibited in response to signaling cause 
derepression of these functions. The convergence of co-regulated 
miRNAs on cancer cell functions raised the possibility that they 
interact to mediate oncogenic control of tumor growth and 
invasion. 

miRNAs regulated by B-Raf/MKK/ERK cooperate in repressing 
invasion 

To test the hypothesis that coordinately regulated miRNAs with 
shared functions might cooperate in augmenting phenotypic 
responses, spheroid invasion was assayed in WM239A cells 
transfected with miRNA mimics that individually repressed 
invasion and proliferation (let-7i, miR-21 1 and miR-22; Figures 3a 
and c), as well as pairwise combinations of these miRNA mimics. 
The combinations of let-7i + miR-22, let-7i + miR-21 1 and let-7i + 
miR-22 + miR-21 1 significantly augmented the effects of each 
miRNA alone by further repressing spheroid cell invasion (Figures 
4a and b). However, no differences were observed with other 
miRNA pairs, indicating that combinatorial repression of spheroid 
invasion was not due to non-specific toxicity (Supplementary 
Figure S5). The results show that miRNA co-regulation by B-Raf/ 
MKK/ERK may augment the effects of individual miRNAs and 
suggest that miRNAs work combinatorially to modulate cancer 
phenotypes. 

miRNAs converge on cancer regulatory genes 
We next considered gene products which might be controlled by 
B-Raf/MKK/ERK-regulated miRNAs, examining cases where a single 
miRNA targets multiple genes and cases where miRNAs converge 
on one gene. TargetScan and miRanda prediction algorithms 23,24 
were used to identify potential targets based on seed sequence 
matches, conservation between species and other criteria. 



However, predicted miRNA targets often include many false 
positives. Therefore, we filtered predicted gene target lists by 
matching them against B-Raf/MKK/ERK-responsive genes 
identified by DNA microarray analysis of the six melanoma cell 
lines treated with U0 1 26. 25 By filtering predicted targets against 
mRNAs responsive to U0126, candidate genes were identified with 
known roles in melanoma proliferation, adhesion/invasion, 
signaling, survival and melanogenesis (Figure 5a). 

The results predicted that gene targets are controlled by 
multiple network miRNAs. Most candidate genes contained 
seed match binding sites for more than one miRNA regulated in 
the same direction, as well as opposite directions, in response to 
B-Raf/MKK/ERK. For example, cyclins D and E were predicted 
targets of miRNAs that were downregulated (let-7, miR-34, -211, 
-26 and -30) as well as upregulated (miR-17/93 and miR-92/25 
families) by B-Raf/MKK/ERK signaling. This suggested that 
while inhibiting miRNAs individually may induce cyclin, the 
effect of many pathway-regulated miRNAs might be more 
complex. We therefore asked whether co-regulated and 
oppositely regulated miRNAs within the B-Raf/MKK/ERK network 
converge on genes. 

To explore this question, we monitored the effect of different 
miRNA mimics on eight proteins, each confirmed a target of B-Raf/ 
MKK/ERK based on responses to two MKK1/2 inhibitors, U0126 
and CI-1040 (Figure 5b). The eight genes had 29 predictions for 
targeting by B-Raf/MKK/ERK miRNAs, of which 22 (-75%) 
accurately predicted protein repression. Thus, filtering against 
cDNA microarray measurements yielded reasonably accurate 
target identification. Proteins were coordinately repressed by 
members of seed-matched miRNA families, such as miR-17/20a/ 
106a/106b/93 and let-7i/7g (Figure 5b), consistent with inhibitory 
mechanisms involving shared 3'UTR seed sequences. Five proteins 
(CCND1, p21, ITGB3, LEF1 and FOSL1) were repressed by two or 
more predicted miRNAs that were downregulated by B-Raf/MKK/ 
ERK, confirming that co-regulated miRNAs converge on proteins. 
Additionally, all six genes containing 3'UTR seed sites correspond- 
ing to oppositely regulated miRNAs were indeed repressed by at 
least one downregulated and one upregulated miRNA. Examina- 
tion of luciferase reporters containing 3'UTR sequences from 
CCND1, CCNE2 and ITGB3 confirmed that each gene was directly 
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Figure 3. miRNAs regulated by B-Raf/MKK/ERK control proliferation and 3D invasion of melanoma cells, (a) WM239A cells transfected with 
1 00 nM miRNA were assayed for cell proliferation using WST-1 . Values 3 days after plating were normalized to values at the time of plating, and 
then normalized to DMRIE-C-treated control cells. NT: negative control (non-targeting) miRNA. Error bars represent s.e.m. (n = 3). (b) Cell 
confluency after miRNA transfection. Phase microscopy images were collected on cells transfected as in (a), after 4 days, (c, d) Spheroid 
growth and invasion following transfection of miRNAs (c) upregulated or (d) downregulated by B-Raf/MKK/ERK. 



targeted by miRNAs that were repressed as well as induced by 
B-Raf/MKK/ERK (Figure 5c). Together, the results show that 
multiple network miRNAs converge on 3'UTRs, suggesting that 
they might interact to control protein repression. 



Combinatorial miRNA expression reveals complex effects on 
protein expression 

Given the evidence for coordinate regulation of multiple proteins, 
we asked how B-Raf/MKK/ERK regulated miRNAs work in 
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Figure 4. Co-expression of miRNAs enhances suppression of cell invasion, (a) Spheroid invasion following transfection of miRNAs individually 
and in combination. Cells were transfected with 50nivi of each miRNA mimic and imaged for cell invasion 3 days after collagen implantation. 
Duplicate spheroids after 3 days in collagen are shown, (b) Cell invasion was quantified by the difference between the outer (R2) and inner 
(R1) spheroid radii (Supplementary Figure S3, panel a), and normalized to control spheroids transfected with negative control (NT: non- 
targeting) miRNA. Asterisks indicate individual miRNAs differing significantly from NT (P-value < 0.003). Double crosses Q) designate miRNA 
combinations that differ significantly from each individual miRNA (P-value < 0.05). Error bars represent s.e.m. (n = 3). 



combination to control protein expression. Four genes were examined in cells transfected with individual miRNA mimics or 
selected for analysis, based on their miRNA binding sites separated combinations of two miRNA mimics, where one, both or neither 
by varying distances (Figure 6a). ITGB3, p21, LEF1 and FOSL1 were miRNA in the combination contained a direct seed match to the 
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Figure 5. miRNAs within the B-Raf/MKK/ERK network cross-regulate cancer regulatory proteins, (a) Cancer genes controlled by B-Raf/MKK/ERK 
were matched to gene targets of miRNAs predicted by either TargetScan (conserved and non-conserved) or miRanda. (b) Western blots show 
repression of proteins corresponding to selected cancer genes, following WM239A transfection with miRNA mimics (50 nM) for 48 h. TBP 
served as the loading control. Underlined miRNAs designate family members with similar seed sequence, which are therefore likely to share 
targets. Asterisk indicates two upper bands of LEF1 which represent longer variants that contain the 3'UTR in the transcript. The lower two 
bands represent variants lacking the 3'UTR. Vertical lines indicate where lanes have been removed from the blots for clarity, (c) Luciferase 
reporter assays validate CCND1, ITGB3 and CCNE2 as direct targets of selected miRNAs. 3'UTR luciferase activity was measured 48 h after 
reporters were co-transfected with 20 nM mimic. Values were normalized first to reporters treated with a B-Raf/MKK/ERK miRNA mimic with no 
binding site in the 3'UTR (NT) and then to cells transfected in parallel with GAPDH-3'UTR reporter. Error bars represent s.e.m. of triplicates. 
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Figure 6. Combinatorial effects on protein repression by miRNAs within the B-Raf/MKK/ERK network, (a) Predicted miRNA binding sites within 
the 3'UTR of four genes, (b) Western blots show protein expression following treatment with miRNA mimics, individually and in combination. 
WM239A cells were transfected with 15 nivi of each miRNA and/or negative control (NT: non-targeting) miRNA (30nM total). qRT-PCR analyses 
of WM239A cells transfected with individual miRNAs show fold increases in expression of ~ 40-200, verifying that the miRNA mimics increase 
the levels of mature miRNA (data not shown), (c) Quantification of western blots in (b), from cells transfected in biological triplicates, 
normalized to TBP and then NT control mimic. Data show mean and s.e.m. (n = 3). Arrows (|) indicate miRNAs predicted to bind 3'UTRs. 
(d) Luciferase reporter assays for miRNA combinations. 3'UTR luciferase activity was measured 48 h after WM239A cells were co-transfected 
with luciferase reporter constructs and 20nivi total miRNA mimics (1 0 r»M each in combination, 1 0 nivi each plus 10nivi non-targeting for 
individual miRNAs). Luciferase activity for each gene was normalized to a control reporter (lacking a 3'UTR) treated in parallel, and then to 
non-targeting control treatment. For (c) and (d), asterisks (*) specify proteins or 3'UTR activity repressed by individual miRNAs vs non-targeting 
control (P-value <0.05). Single crosses (f) indicate miRNA combinations differing from individual targeting miRNAs (P-value <0.07). Double 
crosses (J) indicate miRNA combinations differing from both individual miRNAs (P-value <0.05). 



gene 3'UTR (Figure 6b). In general, individual miRNAs repressed only 
proteins containing at least one 3'UTR seed sequence match 
(Figure 6c, arrows and asterisks) and did not affect proteins without 
seed matches. Exceptions included the lack of FOSL1 repression by 
seed-matched miR-17, and the observed p21 repression by miR-34a 
although no seed matches are present (Figures 6a and c). 

Protein expression analysis after co-transfection with two different 
miRNA mimics revealed varied interactions (Figure 6c; Supplementary 



Figure S6). Generally, co-expressing two 'direct targeting' miRNAs 
(meaning each miRNA has a matched seed sequence and repressed 
protein expression) resulted in augmented protein repression 
compared with individual miRNAs. For example, ITGB3 was further 
repressed by the let-7i + miR-92 combination compared with let-7i or 
miR-92 alone, and LEF-1 was further repressed by miR-1 7 + miR-22 
compared with miR-1 7 or miR-22 alone (Figure 6c, double crosses 'X). 
On the other hand, no further repression of LEF1 or FOSL1 was 
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observed when miR-34 was combined with other direct targeting 
miRNAs (Figure 6c), as if the repression of LEF1 and FOSL1 by miR-34a 
alone was already saturating. 

Combinations of a direct targeting miRNA and a 'non-targeting' 
miRNA (meaning did not repress protein expression) showed 
complex effects. In most cases, protein repression resembled 
the direct targeting miRNA alone (Figure 6c, combinations 
denoted by 'o'). For example, ITGB3 repression by the combination 
of any non-targeting miRNA with either let-7i or miR-92a equaled 
that of let-7i or miR-92a alone. Likewise, p21 repression by the 
combination of any non-targeting miRNA with miR-17 or miR-22 
equaled that of miR-17 or miR-22 alone. Thus, most non-targeting 
miRNAs yielded no interactions. However, we noted exceptions to 
this rule. Although miR-17 individually did not affect FOSL1, it 
augmented repression beyond that seen with miR-22 alone 
when the two miRNAs were combined (Figure 6c). This suggests 
that certain miRNAs show a 'latency' effect, where protein 
repression is observed only when co-expressed with another 
miRNA. Unexpectedly, we also observed interference by a non- 
targeting miRNA toward targeting miRNAs (Figure 6c, single 
crosses T). For example, the non-targeting let-7i derepressed 
LEF1 and FOSL1 proteins when combined with direct-targeting 
miRNAs (miR-17, -22, -34a for LEF1; miR-22, -34a for FOSL1). 
Augmentation and interference effects between miRNAs were 
also observed for 3'UTR activity upon co-transfection of luciferase 
reporters and miRNA combinations (Figure 6d), suggesting that 
these interactions may occur at the level of the 3'UTR. Altogether, 
the results demonstrate that protein expression is influenced by 
complex rules that govern miRNA-specific and gene-specific 
interactions between coordinately regulated miRNAs. 

Knockdown of gene targets phenocopies miRNA effects on cell 
proliferation and invasion 

Finally, we asked whether protein targets regulated by B-Raf/MKK/ 
ERK miRNAs could contribute to miRNA effects on cell prolifera- 
tion or invasion. FOSL1 (miR-34a target) and ITGB3 (let-7i target) 
expression was reduced using siRNA treatment (Figure 7a) and cell 
proliferation and/or invasion was analyzed. Both FOSL1 and ITGB3 
siRNA repressed cell proliferation, similarly to levels seen with miR- 
34a or let-7i (Figure 7b). Additionally, ITGB3 knockdown also 
phenocopied let-7i invasion effects (Figure 3d; Supplementary 
Figure S3) by repressing spheroid invasion to comparable levels 
(Figures 7c and d). The results suggest that complex regulation of 
ITGB3 and FOSL1 by B-Raf/MKK/ERK miRNAs may contribute to the 
widespread influence of the network miRNAs on oncogenic 
melanoma cell behaviors. 



DISCUSSION 

Using a sensitive miRNA array platform, we established a network 
of >20 miRNAs regulated by oncogenic B-Raf signaling in 
melanoma cells. Importantly, the majority of miRNAs controlled 
melanoma cell growth, invasion and gene expression. Interactions 
between co-expressed miRNAs augmented repression of cancer- 
associated proteins and cell invasion from spheroids. We also 
observed unexpected, complex interactions between co- 
expressed miRNAs. Altogether, our findings reveal that molecular 
targets of oncogenic B-Raf include an miRNA network, which 
controls cancer cell properties, and network members in 
combination show different cross-regulatory interactions. 

The miRNAs controlled by B-Raf/MKK/ERK in our screen were 
consistent with those reported in melanoma tumors and cell lines. 
For example, miR-211 and let-7i, which were repressed across 
melanoma cells in our study, were also repressed in melanoma 
tumors compared with melanocytes. 12,15 Likewise, our results 
agreed with reports that miRs-10b, 34a, 34b, 22, 125a and 132 are 
significantly downregulated in primary and metastatic melanoma 
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cell lines. 12 We also observed increased miR-221 and miR-222 
expression, which was previously reported for melanoma 
tumors. 26 Finally, six miR-17 ~92 family and cluster members 
(miR-17-5p, miR-20a, miR-106b, miR-106a, miR-106b and miR-92b) 
were significantly elevated in our melanoma cells, consistent with 
the upregulation of family members in melanoma tissues. 15,27 Our 
results demonstrate that miRNAs previously shown to be 
deregulated in melanomas are controlled by oncogenic B-Raf/ 
MKK/ERK signaling. 

It was intriguing that a large number of miRNAs were co- 
regulated in response to one signaling pathway. Significantly, we 
found that most miRNAs repressed by B-Raf/MKK/ERK inhibited 
cell proliferation and invasion, while miRNAs induced by this 
pathway enhanced these responses. Thus, the directionality of 
miRNA regulation by oncogenic B-Raf/MKK/ERK signaling was 
consistent with promoting cancer phenotypes. This raised the 
question of whether miRNAs coordinately repressed or induced by 
oncogenic signaling are redundant, or whether they interact to 
augment cancer cell behavior. The latter hypothesis was 
corroborated by experiments demonstrating cooperative interac- 
tions between let-7i, miR-22 and miR-211 in suppressing spheroid 
cell invasion. It appears that cell invasiveness in melanoma can be 
augmented by combinatorial interactions between coordinately 
regulated miRNAs within the network. 

Another interesting finding was that coordinately regulated 
miRNAs often converge on the same gene target, which allowed 
us to explore combinatorial miRNA effects on protein expression 
for several genes. In most cases, co-transfecting two miRNAs 
that had 3'UTR seed matches and individually repressed protein 
levels resulted in further repression. This was consistent with 
prior studies demonstrating multiplicative effects on transcript 
abundance when 3'UTR regions contained multiple seed matches 
to the same or different miRNAs. 28,29 We observed novel, 
more complex interactions upon co-transfection with direct- 
targeting and non-targeting miRNAs. The augmented repression 
in the combination of miR-17 + miR-22 revealed a latency effect 
by miR-17 that was gene specific, because unlike FOSL1, miR-17 
targeted both p21 and LEF1 when added alone and in 
combination with miR-22 (Supplementary Figure S6). We also 
observed interference effects, where the ability of one miRNA to 
repress protein was blocked by another. To our knowledge, this is 
the first demonstration that miRNAs without overlapping binding 
sites are able to cancel the effects of others. Conceivably, 
interference between miRNAs could be explained by competition 
for cell uptake or binding the Ago-RISC complex. However, we 
believe this unlikely, because the miRNAs showed interference 
toward some proteins but not others in the same co-transfected 
samples. Altogether, our results reveal unexpected combinatorial 
interactions between components of an oncogene-regulated 
miRNA network. The interactions are specific for miRNAs and 
proteins; thus, miRNA regulatory mechanisms are complex 
and rules for combinatorial regulation cannot be assumed. 

Our findings show that B-Raf/MKK/ERK may have far-reaching 
effects in regulating cancer cell behavior through the control 
of miRNA expression. Each miRNA within the network targets 
many genes which could contribute to the cancer phenotype 
(Figure 8). Importantly, our study uncovers new combinatorial 
interactions between miRNAs. Coordinate regulation of an miRNA 
network by a signaling pathway may lead to unpredictable effects 
on proteins when combinatorial effects are considered, and 
further exploration of the rules for such interactions are needed. 
Finally, while our study focuses on only one signaling pathway, 
most cancers (including melanoma) harbor many deregulated 
pathways. To the extent that different signaling pathways each 
control different miRNA networks, an important question for 
future investigation is whether these networks interact combina- 
torially to relay the effects of multiple pathways on the cell 
proteome. 
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Figure 7. Knockdown of FOSL1 and ITGB3 phenocopy their targeting miRNAs. (a) Western blots show repression of FOSL1 and ITGB3 
expression 48 h after WM239A cells were transfected with 50nivi siRNA. Two individual duplexes (dup 1, dup 2) were used for each gene. TBP 
serves as loading control, (b) Cell count proliferation assays for WM239A cells transfected with 50nivi miRNA mimics or siRNAs. In all, 100 000 
cells were plated in triplicate on 6-well plates and counted 5 days after plating. Cell counts were normalized to cell concentrations obtained 
using the WST-1 assay at the time of plating, (c) Spheroid invasion was analyzed for WM239A cells treated with 50nivi control siRNA or ITGB3 
siRNA 48 h before spheroid formation, (d) Spheroid invasion was quantified by the difference between the outer (R2) and inner (R1) spheroid 
radii, and normalized to control spheroids transfected with control, non-targeting siRNA. Error bars represent s.e.m. for three spheroids. 
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Figure 8. Model illustrating part of the miRNA network regulated by 
B-Raf/MKK/ERK in melanoma cells. B-Raf/MKK/ERK induces or 
represses sets of miRNAs, each which target melanoma-associated 
oncogenes or tumor suppressors. Multiple co-regulated miRNAs 
may converge on and cross-regulate each gene, allowing oncogenic 
B-Raf to exert variable regulatory effects on melanoma genesis. 



MATERIALS AND METHODS 

Cell culture 

Primary melanocyte cell lines NHEM693 and NHEM2493 (Cambrex Bio 
Sciences, Walkersville, MD, USA) were cultured in Medium154 plus Human 
Melanocyte Growth Supplement (Cascade Biologies, Portland, OR, USA). 
Melanoma cell lines were a gift from Meenhard Herlyn, Wistar Institute, 30 
and were grown in RPMI (Invitrogen, Carlsbad, CA, USA) supplemented 



with 10% FBS (Gemini, West Sacramento, CA, USA) and penicillin- 
streptomycin (100U/ml). All cell lines were maintained at 37 °C and 

5% C0 2 . 



miRNA microarray profiling 

Melanocyte and melanoma cell lines were plated on 6 cm plates to achieve 
~60% confluency after 48 h, at which time cells were harvested in Buffer 
RLT Plus (Qiagen, Germantown, MD, USA) and stored at -80°C For 
melanoma cell drug treatments, 1 0 ja/vi U01 26 (Promega, Madison, Wl, USA) 
or 3(im CI-1040 (Pfizer, Inc., New York, NY, USA) in 0.1% (v/v) dimethyl 
sulfoxide or 0.1% dimethyl sulfoxide carrier was added 24 h before 
harvesting. Biological replicates were generated from different cell 
passages, each ^2 passages removed. Total RNA was purified using the 
RNeasy Plus kit (Qiagen) according to manufacturer's Supplementary 
Protocol 1. MicroRNAs were quantified using the Thermo Scientific 
Dharmacon microRNA expression profiling platform, which uses two-color 
high-density eight-plexed slides containing probes to 470 mature human 
miRNAs. Total RNA (500 ng) was combined with spike-in control to monitor 
labeling efficiency, then fluorescence dye-labeled with pCp-Cy5 using T4 
RNA ligase (Promega). A control oligonucleotide was added to measure 
recovery, and the sample desalted through a P6 Biospin column (Bio-Rad, 
Hercules, CA, USA). The miRIDIAN Synthetic miRNA Reference Probe Set 
(Thermo Fisher Scientific), 2 x HI-RPM hybridization buffer (Agilent, Palo 
Alto, CA, USA) and 25 x fragmentation buffer were added to each sample 
and hybridized to eight 15 K custom miRNA microarrays (miRBase 9.0, 
probes in triplicate) (Agilent). Hybridizations were performed in a rotating 
incubator at 1 2 r.p.m. and 53 °C for 20 h. Slides were washed and scanned 
on an Agilent Microarray Scanner (Model G2505B), and data were 
extracted using Agilent Feature Extraction 9.5.1 software. Triplicate probe 
signals were averaged and adjusted for background, and the Cy5 (sample) 
channel was normalized to Cy3 (synthetic reference) for each probe. For 
each array, the normalized Cy5 signals were scaled against spike-in 
controls to adjust for global labeling bias. 'Is-Present' P-values were 
calculated for signals in the Cy5 channel using a set of background (non- 
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specific hybridizing) probes and the Gaussian error function. A set of 264 
high-confidence probes, defined as those with P<0.05 in 20 of 40 samples, 
was used as the basis feature set for inter-array normalization, using a 
mean-centering scaling. The normalized signals were transformed to log qo 
and analyzed using Microsoft Excel. 

Real-time quantitative RT-PCR 

Cell lysates were prepared as described for microarray analysis, and real- 
time RT-PCR by absolute quantification was performed by contract 
(Asuragen, Inc., Austin, TX, USA). Briefly, total RNA was isolated using a 
column-based method and quality was assessed with an Agilent 2100 
Bioanalyzer (Agilent). In all, 1 ng of total RNA was reverse transcribed using 
ABI Taqman miRNA assays (Applied Biosystems, Austin, TX, USA) and 3 uJ of 
the product was used in subsequent Taqman PCR amplification, performed 
in triplicate on a validated real-time thermocycler (ABI 7900HT). C t values 
were compared with standard curves generated with synthetic miRNAs to 
obtain absolute copy number per ng RNA, and values were normalized 
using stably expressed RNU44 small nucleolar RNA. 

siRNA, miRNA oligonucleotide transfections 

Synthetic oligonucleotide siGENOME siRNAs and miRIDIAN miRNA Mimics 
and Hairpin Inhibitors (Thermo Fisher Scientific) were reconstituted in 1 x 
sterile siRNA buffer (60 mM KCI, 6mivi HEPES-pH 7.5, 0.2mivi MgCI 2 ) to a 
stock concentration of 20ujvi. Cells were transfected in suspension using 
DMRIE-C (Invitrogen) and incubated for 5 h before replacing the media 
with RPMI + 10% FBS. 

WST-1 cell proliferation assay 

Cells were transfected with miRNA mimic or DMRIE-C lipid only. After 24 h, 
5000 cells were plated in triplicate on duplicate 96-well plates. Three hours 
before reading, a standard curve of 500-50000 cells was added. Plates 
were centrifuged at lOOOr.p.m. for 5 min before replacing the media with 
100uJ media containing WST-1 reagent at 60uJ/ml (Clontech, Mountain 
View, CA, USA). Plates were incubated at 37 °C for 1 h and absorbances 
quantified at 440 nm. 

Spheroid invasion assays 

Spheroids were grown in collagen matrix as described. 31 Twenty-four 
hours after miRNA transfection, cells were plated in triplicate on 1.5% 
Noble agar coated 96-well plates, then incubated for 1-2 days. Spheroids 
were collected and centrifuged, resuspended in 2.5 mg/ml PureCol bovine 
collagen-l (INAMED, Leimuiden, Netherlands) containing RPMI + 10% FBS, 
and overlaid on a pre-solidified layer of the same collagen solution in 
96-well plates. The collagen/spheroid suspension was solidified for 1 h, 
then overlaid with RPMI + 10% FBS. Microscopy images of spheroid growth 
and invasion were quantified using Image J software (NIH, Bethesda, 
MD, USA). 

Western blotting 

Cells were lysed in 20mivi Tris pH 8.0, 2% SDS, 10mivi dithiothreitol and 
boiled for 10 min. Total protein was determined using the DC Protein Assay 
(Bio-Rad). Membranes were blocked with 5% (w/v) non-fat dry milk and 
washed in phosphate- buffered saline + 0.1% (v/v) Tween-20. Primary 
antibodies were diluted in 5% BSA + phosphate- buffered saline/0.1% 
Tween-20 and recognized ERK1/2, LEF-1, CCNE2, p21, FOSL1 (Cell Signaling 
Technology, Danvers, MA, USA), MLANA, ITGB3, DUSP4, TBP (Santa Cruz 
Biotechnology, Santa Cruz, CA, USA), phospho-ERK1/2 (Sigma, St Louis, MO, 
USA) and CCND1 (Upstate, Charlottesville, VA, USA). Horseradish perox- 
idase-coupled secondary antibodies (Jackson ImmunoResearch, West 
Grove, PA, USA) were diluted 1:5000 in phosphate-buffered saline + 0.1% 
Tween-20 with 5% non-fat milk and visualized using enhanced chemilu- 
minescence (Perkin-Elmer, Waltham, MA, USA). 

3'UTR luciferase reporter assays 

Melanoma cells were co-transfected with miRNA mimics and 3'UTR 
luciferase reporters (Switchgear Genomics, Menlo Park, CA, USA, individual 
miRNA assays; or GeneCopoeia, Rockville, MD, USA, combinatorial miRNA 
assays) using DharmaFECT Duo (Thermo Fisher Scientific). In all, 20000 
cells/well were plated onto opaque 96-well plates (Becton Dickinson, 
Franklin Lakes, NJ, USA), incubated for 24 h, and transfected with lOOng 
(SwitchGear Genomics) or 50 ng (GeneCopoeia) 3'UTR luciferase reporter 
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plus 5-20 nivi mimic. Luminescence was analyzed after 24 h by replacing 
media with 75 |il Steady-Glo luciferase reagent (Promega), and measuring 
luminescence after 30 min (Beckman Coulter DTX 880, Fullerton, CA, USA). 

CONFLICT OF INTEREST 

Dr Emily Anderson, Maren Mayer and Kevin Sullivan are employed by and own stock 
in Dharmacon, Inc., a wholly owned subsidiary of Thermo Fisher Scientific, Inc. 
Remaining authors declare no conflict of interest. 

ACKNOWLEDGEMENTS 

We are indebted to Meenhard Herlyn for melanoma cell lines used throughout these 
studies, and to Judith Leopold-Sebolt for a generous gift of CI-1040. We thank 
SwitchGear Genomics for providing 3'UTR constructs. Finally, we thank Thermo Fisher 
Scientific investigators Amanda Birmingham, Annaleen Vermeulen and Barbara 
Robertson for many helpful discussions. This study was financially supported by NIH 
R01-CA1 18972 (NGA) and NIH T32-GM008759 (KLC). 



REFERENCES 

1 Davies H, Bignell GR, Cox C, Stephens P, Edkins S, Clegg S et al. Mutations of the 
BRAF gene in human cancer. Nature 2002; 417: 949-954. 

2 Edlundh-Rose E, Egyhazi S, Omholt K, Mansson-Brahme E, Platz A, Hansson J et al. 
NRAS and BRAF mutations in melanoma tumours in relation to clinical char- 
acteristic: a study based on mutation screening by pyrosequencing. Melanoma 
Res 2006; 16: 471-478. 

3 Koo HM, VanBrocklin M, McWilliams MJ, Leppla SH, Duesbery NS, Woude GF. 
Apoptosis and melanogenesis in human melanoma cells induced by anthrax 
lethal factor inactivation of mitogen activated protein kinase kinase. Proc Natl 
Acad Sci USA 2002; 99: 3052-3057. 

4 Hingorani SR, Jacobetz MA, Robertson GP, Herlyn M, Tuveson DA. Suppression of 
BRAF(V599E) in human melanoma abrogates transformation. Cancer Res 2003; 63: 
5198-5202. 

5 Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 2004; 
116: 281-297. 

6 Lee YS, Dutta A. MicroRNAs in cancer. Annu Rev Pathol 2009; 4: 199-227. 

7 He H, Jazdzewski K, Li W, Liyanarachchi S, Nagy R, Volinia S et al. The role of 
microRNA genes in papillary thyroid carcinoma. Proc Natl Acad Sci USA 2005; 102: 
19075-19080. 

8 Johnson CD, Esquela-Kerscher A, Stefani G, Byrom M, Kelnar K, Ovcharenko D 
et al. The let-7 microRNA represses cell proliferation pathways in human cells. 
Cancer Res 2007; 67: 7713-7722. 

9 Medina PP, Nolde M, Slack FJ. OncomiR addiction in an in vivo model of micro- 
RNA-21 -induced pre-B-cell lymphoma. Nature 2010; 467: 86-90. 

10 Chang TC, Yu D, Lee YS, Wentzel EA, Arking DE, West KM et al. Widespread 
microRNA repression by Myc contributes to tumorigenesis. Nat Genet 2008; 40: 
43-50. 

1 1 Schulte JH, Horn S, Otto T, Samans B, Heukamp LC, Eilers UC et al. MYCN regulates 
oncogenic MicroRNAs in neuroblastoma. Int J Cancer 2008; 122: 699-704. 

12 Caramuta S, Egyhazi S, Rodolfo M, Witten D, Hansson J, Larsson C et al. MicroRNA 
expression profiles associated with mutational status and survival in malignant 
melanoma. J Invest Dermatol 2010; 130: 2062-2070. 

13 Chan E, Patel R, Nallur S, Ratner E, Bacchiocchi A, Hoyt K et al. MicroRNA 
signatures differentiate melanoma subtypes. Cell Cycle 2011; 11: 1845-1852. 

14 Mueller DW, Rehli M, Bosserhoff AK. miRNA expression profiling in melanocytes 
and melanoma cell lines reveals miRNAs associated with formation and pro- 
gression of malignant melanoma. J Invest Dermatol 2009; 129: 1740-1751. 

15 Chen J, Feilotter HE, Pare GC, Zhang X, Pemberton JG, Garady C et al. MicroRNA- 
193b represses cell proliferation and regulates cyclin D1 in melanoma. Am J 
Pathol 2010; 176: 2520-2529. 

16 Philippidou D, Schmitt M, Moser D, Margue C, Nazarov PV, Muller A et al. 
Signatures of microRNAs and selected microRNA target genes in human 
melanoma. Cancer Res 2010; 70: 4163-4173. 

17 Satyamoorthy K, Li G, Gerrero MR, Brose MS, Volpe P, Weber BL et al. Constitutive 
mitogen-activated protein kinase activation in melanoma is mediated by both 
BRAF mutations and autocrine growth factor stimulation. Cancer Res 2003; 63: 
756-759. 

18 Mills CN, Joshi SS, Niles RM. Expression and function of hypoxia inducible factor-1 
alpha in human melanoma under non-hypoxic conditions. Mol Cancer 2009; 
8: 104. 

19 Paroo Z, Ye X, Chen S, Liu Q. Phosphorylation of the human microRNA-generating 
complex mediates MAPK/Erk signaling. Cell 2009; 139: 112-122. 

20 Mendell JT. miRiad roles for the miR-17-92 cluster in development and disease. 
Cell 2008; 133: 217-222. 



© 2013 Macmillan Publishers Limited 



Oncogene (2013) 1959-1970 



B-Raf controls a network of microRNAs in melanoma 

KL Couts et al 



1970 

21 Roush S, Slack FJ. The let-7 family of microRNAs. Trends Cell Biol 2008; 18: 505-516. 

22 Hermeking H. The miR-34 family in cancer and apoptosis. Cell Death Differ 2010; 
17: 193-199. 

23 Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adeno- 
sines, indicates that thousands of human genes are microRNA targets. Cell 2005; 
120: 15-20. 

24 Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS. Human MicroRNAtargets. PLoS 
Biol 2004; 2: e363. 

25 Argast GM, Croy CH, Couts KL, Zhang Z, Litman E, Chan DC et al. Plexin B1 is 
repressed by oncogenic B-Raf signaling and functions as a tumor suppressor in 
melanoma cells. Oncogene 2009; 28: 2697-2709. 

26 Felicetti F, Errico MC, Bottero L, Segnalini P, Stoppacciaro A, Biffoni M et al. The 
promyelocytic leukemia zinc finger-microRNA-221/-222 pathway controls 
melanoma progression through multiple oncogenic mechanisms. Cancer Res 
2008; 68: 2745-2754. 

27 Igoucheva O, Alexeev V. MicroRNA-dependent regulation of cKit in cutaneous 
melanoma. Biochem Biophys Res Commun 2009; 379: 790-794. 

28 Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA 
targeting specificity in mammals: determinants beyond seed pairing. Mol Cell 
2007; 27: 91-105. 

29 Saetrom P, Heale BS, Snove Jr O, Aagaard L, Alluin J, Rossi JJ. Distance constraints 
between microRNA target sites dictate efficacy and cooperativity. Nucleic Acids 
Res 2007; 35: 2333-2342. 

30 Smalley KS, Haass NK, Brafford PA, Lioni M, Flaherty KT, Herlyn M. Multiple 
signaling pathways must be targeted to overcome drug resistance in cell lines 
derived from melanoma metastases. Mol Cancer Ther 2006; 5: 1136-1144. 

31 Hsu MY, Elder DE, Herlyn M. Melanoma: the Wistar melanoma (WM) cell lines. 
Human Cell Culture 1999; 1: 259-274. 

32 Yu F, Yao H, Zhu P, Zhang X, Pan Q, Gong C et al. let-7 regulates self renewal and 
tumorigenicity of breast cancer cells. Cell 2007; 131: 1109-1123. 



33 He L, He X, Lowe SW, Hannon GJ. microRNAs join the p53 network-another piece 
in the tumour-suppression puzzle. Nat Rev Cancer 2007; 7: 819-822. 

34 Xiong J, Yu D, Wei N, Fu H, Cai T, Huang Y et al. An estrogen receptor alpha 
suppressor, microRNA-22, is downregulated in estrogen receptor alpha-positive 
human breast cancer cell lines and clinical samples. FEBS J 201 0; 277: 1 684-1 694. 

35 Huse JT, Brennan C, Hambardzumyan D, Wee B, Pena J, Rouhanifard SH et al. The 
PTEN-regulating microRNA miR-26a is amplified in high-grade glioma and 
facilitates gliomagenesis in vivo. Genes Dev 2009; 23: 1327-1337. 

36 Braun J, Hoang-Vu C, Dralle H, Huttelmaier S. Downregulation of microRNAs 
directs the EMTand invasive potential of anaplastic thyroid carcinomas. Oncogene 
2010; 29: 4237-4244. 

37 Landais S, Landry S, Legault P, Rassart E. Oncogenic potential of the miR-1 06-363 
cluster and its implication in human T-cell leukemia. Cancer Res 2007; 67: 
5699-5707. 

38 Ambs S, Prueitt RL, Yi M, Hudson RS, Howe TM, Petrocca F et al. Genomic profiling 
of microRNA and messenger RNA reveals deregulated microRNA expression in 
prostate cancer. Cancer Res 2008; 68: 6162-6170. 

39 Li Y, Tan W, Neo TW, Aung MO, Wasser S, Lim SG et al. Role of the miR-1 06b-25 
microRNA cluster in hepatocellular carcinoma. Cancer Sci 2009; 100: 1234-1242. 

40 Nagaraja AK, Creighton CJ, Yu Z, Zhu H, Gunaratne PH, Reid JG et al. A link 
between mir-100 and FRAP1/mTOR in clear cell ovarian cancer. Mol Endocrinol 
2010; 24: 447-463. 

41 Papadopoulos GL, Reczko M, Simossis VA, Sethupathy P, Hatzigeorgiou AG. The 
database of experimentally supported targets: a functional update of TarBase. 
Nucleic Acids Res 2009; 37: D1 55-D1 58. 



0 | This work is licensed under the Creative Commons Attribution- 
Wi i i i i i i ii H i ira mjiJ NonCommercial-Share Alike 3.0 Unported License. To view a copy of 
this license, visit http://creativecommons.Org/licenses/by-nc-sa/3.0/ 



Supplementary Information accompanies the paper on the Oncogene website (http://www.nature.com/onc) 



Oncogene (2013) 1959-1970 



© 2013 Macmillan Publishers Limited 



